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£^ ■ The pion electromagnetic form factor is investigated by using the dispersion relation with 

the superconvergence condition, which was proposed to synthesize the vector meson dominance 
model and QCD. The absorptive part is given as an addition of the resonance and the QCD 
terms, the latter being parametrized so as to realize the prediction of perturbative QCD. 

The experimental data are analyzed and the parameters for the QCD part, the resonance 
part as well, are determined. Existing experiments for the form factor are reproduced very well 
both for the space-like and time-like momenta. 



<D . 1 Introduction 



PACS 11.55.Fv - Dispersion relations 

PACS 12. 38. Aw - General properties of QCD (dynamics confinement, etc.) 
^ ' PACS 13.40Gp - Electromagnetic form factor. 

p 

According to the perturbative QCD hadron electromagnetic form factors decrease more rapidly than 
an inverse polynomial for large squared momentum transfer [Q, || . The boson electromagnetic form 
H | factor, F(t), decreases as 

F(t)^~327r 2 f 2 B //3 Q tln\t\ (1) 

for \t\ — > oo, where fs is the boson decay constant, t is the squared momentum transfer and (3q is the 
lowest order of the (3 function in the renormalization group, namely /?o = H — 2n//3 with rif being 
the number of flavor. Experimentally, the charge and magnetic form factors of nucleons decrease 
more rapidly for large t than the well-known dipole formula as is predicted by the perturbative QCD. 

The QCD property for the form factor is realized by using the dispersion relation in which the 
absorptive part, Im.F(i), decreases as some power of lnt and satisfies the superconvergence constraint 
H ||. We break up the absorptive part as a simple summation of the resonance and the QCD part. 
The former is given as addition of the Breit-Wigner formula with possible mixing among resonances 
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and the latter as a power series in the running coupling constant which is analytically extrapolated 
to the time-like momentum. We impose the superconvergence condition on lmF(t). For the form 
factor the low momentum part is rather easy to deal with by the dispersion relation. It is, therefore, 
possible to treat whole range, from the low to the high momentum regions, systematically; the vector 
meson dominance model (VMD) is realized for low momentum and the QCD condition is satisfied 
for large \t\. The dispersion relation works as an interpolation. 

In this paper we improve the result of Ref. || , by investigating correlations among QCD parame- 
ters. We search much wider range of parameters than in Ref. Q. We show that QCD part of lmF(t) 
can be determined from the experiments on the pion electromagnetic form factor although the data 
are restricted to low momentum. At least the second order in the running coupling constant is nec- 
essary to fit the experimental data. The QCD scale parameter A is determined to be in the range A 
= 0.3 ~ 2 GeV. The calculated form factor agrees with the existing experimental results very well 
both for the time-like and space-like momenta. 

2 Superconvergence dispersion relation and QCD 

As the superconvergent dispersion relation has already been discussed 0-0, we summarize our 
prescription to make the paper self-contained. To realize the asymptotic form of QCD we make use 
of the following property of the Hilbert transformation. Let F(t) be given as 

F(t) = £dt'p^- t , (2) 

where sq > is a constant and the function p{t') decreases as 

p(t') -> a/{ln(i')} 7 fori'^oo, 
with a and 7 being constants, then 

F(t)= r dt'-^l -> 7 — \ m . fori ^00. 

J so t-t> ( 7 -l){ln(|t|)}7-i 

Here 7 > 1. In the case that the spectral function p(t') — » a/i'ln 7 (t') for t' — ► 00 and satisfies the 

/>OC 

superconvergence condition / p(t')d£ — 0, F(t) given by (Q) decreases asymptotically as F(t) — » 

J so 

a/( 7 - l)t ln 7_1 (|t|) for t -> 00. 

2.1 Vector meson dominance model and QCD 

For the pion electromagnetic form factor, F(t), we assume the unsubtracted dispersion relation with 
respect to the squared momentum transfer t. 

If 00 , ,lmF(t') 
F (t) = ~ / dt'—Y, (3) 

7T Jso 1 ~ 1 
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where sq is the threshold. To satisfy the asymptotic condition of perturbative QCD (Q), the absorp- 
tive part should satisfy the following conditions: 

ImF(t) -» (f -> oo) (4) 

and the superconvergence condition 

/•OO 

dt'Imf((') = 0. (5) 



/ 

J Sri 



' s 

We take ImF so as to realize the result of the vector meson dominance model for low momentum 
region and the QCD prediction asymptotically; we express it as the summation of the Breit-Wigner 
formula with mixing among resonances |l0| and the QCD term. 

lmF(t) = lmF BW (t) + knF QGD (t), (6) 

where the Brcit-Wigner term lmF BW (t) is given as 

(M?-t) 2 + 7 2 

- ^jk{M]-t)+af k l3 oj k {Ml -t) + a? k <y k 

Here the suffixes denote numbering of resonances and a^ k and ol 1 ^ are 

R ljlk {M]-Ml) 



lmF BW (t) = ~t 

i 

^ Cjfe l (M|-t)2+ 7 ? ( Mfc 2 -<) 2 +7 fc 2 J' lj 



ai = - 



(M?-M fe 2 )2 + ( 7j .- 7 ^ 
7i7fe(7j 2 ~7fc) 



(M?-M fe 2 )2 + ( 7j .- 7fe )2 



(8) 
(9) 



Im_F < 5 C£, (t) is given as a power series in the running coupling constant multiplied by a function hit), 
which is introduced to assure the threshold behavior and the convergence of the superconvergence 
integral (||). In fact, 

lmF^ CD (t) = Im [J2 cf° D {a R (t)y] h(t), (10) 
j>i 

Here an is the running coupling constant in QCD analytically regularized || |8| ||] ■ We shall discuss 
on aR in the next subsection. For the function h(t) we assume the following formula 0, motivated 
by Kroll et al.@: 

~t-t ~ 



h(t) = 



3/2 

fc2 



t + V 



(11) 



_t + h_ 

with tj (i = 0, 1, 2) being parameters. It must be noticed that the result is not sensitive to the form 
of the function h(t). 



3 



For the form factor we impose the normalization condition at t = 0, 

F(0) = 1. (12) 

We have two conditions on F(t), the superconvergence and the normalization condition (||) and (p"2|), 
respectively. 



2.2 Running coupling constant for the time-like momentum 

It is known that the running coupling constant as has an unphysical pole for the space-like momen- 
tum if it is calculated by the perturbative QCD. The removal of the singularity was proposed by 
Dokshitzer et al.|^|, ^| via analytic regularization. The regularized running coupling constant, being 
denoted as an, is given by the spectral representation 



1 

ur. = — 

7T 



dt 



v -t 



(13) 



Here the function a is given in terms of the as, which is obtained for the space- like momentum via 
the renormalization group. We have 

1 



= T .[a S {e-™t) - a s {e™t)). 
2,1 

We use the three loop approximation for as, which is expressed by the formula 

4-7T 



(14) 



Po 
+a 2 



ln(Q 2 /A 2 )+ ai ln{ln(g 2 /A 2 )} 
ln{ln(Q 2 /A 2 )} , 
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ln(Q 2 /A 2 ) ln{ln(g 2 /A 2 )} 



where 



Hl = 2/3 1 /f3l a 2 = 4/3 2 //? 4 , a 3 = (4/3 2 //3 4 )(l - /3 /3 2 /8/3 2 ), 
and = 0,1,2) are the /3 functions in QCD with (3q = 11 — 2rty/3, /3i = 51 



(15) 
(16) 

19n//3 and 

/?2 = 2857 — 5033n//9 + 325n 2 /27, nf being the number of flavor. Expanding (15) in terms of 
ln(Q 2 /A 2 ) and ln{ln(Q 2 /A 2 )}, we get the usual expression for the running coupling constant |Q. We 
perform the analytic continuation of the effective coupling constant as to the time-like momentum 
by the replacement Q 2 — > e~ l7! t and decompose to the real and imaginary parts as — Re[as(^)] + 
ilm[as(t)]. Substituting Im[as(t)] in (p^), we obtain the running coupling constant for the time-like 
momentum. We write it by the same notation as an. We showed in Ref.|^] that an given by (13) 
approximately coincides with that which is given by ( |l5| ) with the ghost pole subtracted, namely, 
for the space-like momentum 



a R (Q') « a s (Q 2 ) - A*/(Q 2 - Q* 2 ) 



(17) 



where Q* is the pole of as{Q ) and A* is the residue at the pole. Writing Q* /A = e u , we have 



0.7910487 for nf = 2 with the three loop approximation and 

In it* a 2 



A* = 4nA 2 e u 



/{/3o(l 



ai 

u* 



a 2 - 
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where ai are defined by (ITq) . To extend to time- like momentum we replace Q 2 



r t as is men- 



tioned above. The difference between an given by the spectral representation and the approximate 
one is less than 0.4 % for \t\ £ 3 GeV 2 . The approximation becomes better as \t\ becomes larger 
§. We shall use the approximate one for the regularized running coupling constant hereafter; as we 
have multiplied the function h(t) to define lmF^ CD , the contribution from the low momentum part 
is considerably suppressed so that the approximation turns out to be very good. 



For the time-like momentum an, being obtained by the replacement Q 2 



as follows: 



Re[a R (t)} 
Im[a R (t)} 



47TU 

Pi 

Any 
PoD(t) • 



.4* 



T t in (|17|), is given 

(18) 
(19) 



where 



ln{t/A 2 



a i 



ln{ln(t/A 2 )} 



a 2 



ln 2 (t/A 2 ) +7T 2 

a 3 ln(VA 2 ) 
ln 2 (t/A 2 ) +7T 2 

7r + a\6 — 



{iln(t/A 2 )ln{ln(t/A 2 )} + ^} 



«2 



ln^i/A 2 



— {|ln{ln(t/A 2 )}-01n(t/A 2 )} 



7ra 3 



ln 2 (t/A 2 ) 



2 ' 

r- 



and 



Here 



D = u 2 + v 2 . 



9 = tan 



lln(t/A 2 )/' 

The spectral function a is obtained by taking imaginary part of as, that is, 

cr(t) = 4irv/f3 a D, 

where v and D are defined by (|2l|) and (E2J) , respectively. 



(20) 

(21) 

(22) 
(23) 

(24) 



3 Numerical results 

We analyzed the experimental data of the pion electromagnetic form factor for the space-like and 
the time-like momenta. The parameters appearing in our formulas are determined so as to minimize 
the x 2 , addition of the chi square for the time-like and space-like momenta. We fix the parameters 
which appear in the function hit) at to — A 2 and t\ = t% =16 GeV 2 as the result is not sensitive to 
these values JtJ. The number of flavor is taken as nj = 2. 
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For large t, ImF given by (10), approaches to lmF(t) — > cf CD Ima]i(t)t2 /t and Imajj given by 
(|9|) to 47r 2 //3 ln 2 (t/A). Therefore, comparing @ and (|o|), we have 

cf CD - -87r/|/t 2 , (25) 

which is fixed at c{ CD = —0.02683 by taking = 16 GeV and the pion decay constant f v — 0.1307 
GeV. 

To realize the data for the time-like momentum we take the masses and widths of vector bosons 
as parameters, which are determined as follows: m p = 0.760 GeV, T p = 0.138 GeV; m p > — 1.40 GeV, 
I> = 0.33 GeV; m p „ = 1.73 GeV, T p „ = 0.24 GeV. We investigated the effect of the unconfirmed 
vector boson whose mass and width are taken as m p "i =2.11 GeV and Y p m = 0.368 GeV. The mass 
and width of 10 meson are kept at the experimental values fll2| . We leave out the <j> meson in this 
calculation as the contribution is estimated to be very small. In the chi square analysis the following 
experimental data are used: Amendolia et al.jl4| and Bebek et al.|l5j for the space- like momentum 
and Barkov et al.[^3| for the time-like momentum. The number of data points are 101 for the space- 
like and 157 for the time-like regions; totally 258 data are used in the analysis. Degrees of freedom 
(DOF) is 241 when p'" is omitted. If p'" is taken into account, DOF = 238. 

First, we investigate the QCD scale parameter A. We restrict to the second order approximation 
for the QCD running coupling constant in ImFQ CD , namely, c$ D — 0. We illustrate in Fig.l the 
contour lines for \ 2 = 400, 500 and 600 to give correlation between the pion decay constant /„-, being 
taken as a parameter, and D ■ Here p'" is left out. To draw the contour lines, / ff and c^ CD is kept 
fixed and the other parameters are determined so as to make % 2 minimum. The QCD scale parameter 
is fixed at A = 0.8 GeV. % 2 is then given as a function of f v and c^ CD ■ The horizontal lines in Fig.l 
denote the experimental value for the pion decay constant f„ =0.1307 ± 0.00037 GeV. In the same 
way, we examined correlation between A and c^ CD . In Fig. 2 we show the contour lines which give 
correlation between A, being taken as a parameter, and c^ CD , where the pion decay constant is fixed 
at the experimental value, / w = 0.1307 GeV. Taking \ 2 < 500 in Fig.l and the experimental value 
for /tt, we obtain D — —0.87 ~ —0.49. The one-loop approximation for the QCD part, namely 

° = 0, leads to poor result. The QCD scale parameter is determined from the correlation of 
c% CD and A which is illustrated in Fig.2. The condition x 2 < 500 leads to 0.3 < A < 2 GeV, 
where we take the value of c^ CD determined above, namely, c® CD = —0.87 ~ —0.49. The best fit 
is obtained for A « 0.8 GeV. 
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Fig.l. Contour lines for the parameters c^ CD and f n . \ 2 is obtained as a function of f-n and c^ CD with 
the other parameters determined so as to make x 2 minimum. A is fixed at 0.8 GeV. Solid, dotted and dashed 
lines denote the contours with \ 2 = 400, 500 and 600, respectively. Data for the space-like momentum [ |l4| , 
]15f and the time-like momentum Ml are used. 




Fig. 2. Contour lines for c% and A. The pion decay constant is fixed at the experimental value. See the 
caption of Fig.l. 
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The parameters which minimize \ 2 are summarized in Table I, where A is fixed at A = 0.8 GeV 
and the pion decay constant at the experimental value for the cases with and without p'" . 





without p"'(2.11) 


with p"'{2.n) 


Cp 


2.723 


2.701 




-0.0030 


-0.0030 


c p , 


-0.2701 


-0.0924 


Cp" 


-0.1956 


-0.2830 


Cp,n 




-0.0260 


Cuj'-p 


1.069 


1.068 


Cp-pi 


1.832 


-0.089 


Cp-p" 


0.562 


-0.136 


Cp'-p" 


-51.92 


-49.99 


QCD f 
c l 


-0.02683 


-0.02683 


QCD 
c 2 


-0.645 


-0.573 


x 2 


379.1 


355.1 



f cf CD is calculated by ©. 
Table I Residues at the resonance poles, mixing parameters among resonances and the QCD parameters 
cf CD (i = 1, 2). We take A = 0.8 GeV and cf CD =0. The parameters in the function h(i) are to = A 2 , ti = 
ta = 16 GeV 2 . Two cases are investigated without the unconfirmed resonance p'" and with p'". The value 
of x 2 is the summation of that for the space-like and time-like momenta. In Figs. 1-5 the same parameters 
are used. 

Let us compare our calculated results with the experimental data of the pion electromagnetic 
form factor. In Fig. 3 we give our results for ^(t)! 2 for the space-like momentum. The experimental 
data are taken from Amendolia et al.[Q. In Fig. 4 \t \ F(t) is compared with the experiments. Our 
result agrees with the recent experiment by J.Volmer et al.[|l6). In the figure the data by Bebek et 
al.|l5|] are included. 
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'o 0.1 0.2 0.3 

-t (GeV 2 ) 



Fig.3. Electromagnetic form factor |_F(t)| 2 for the space-like momentum. Low momentum region. Data are 
taken from Amendolia et al. |p^| . Two cases with and without p" almost coincide. The parameters given in 
Table I are used. 



t 1 1 1 1 1 1 1 1 r 




-t (GeV 2 ) 

Fig. 4. t F(t) for the space- like momentum. The data denoted as x are taken from Bebek et al.Jl^], the open 
circles and closed circles are taken from Volmer at al. fllq] . Solid line stands for the calculation without p'" 
and the dashed one with p" . The parameters given in Table I are used. 



The solid lines stand for the calculation without assuming p'" and the dashed ones with p 1 




Fig. 5 \F(t)\ 2 for the time-like momentum. Data are taken from Barkov et al.|Q|. The solid lines stand for 
the calculation without p" and the dashed ones with p'" . (a) t = 0.3 - 0.75 GeV 2 , (b) t = 0.75 - 0.85 
GeV 2 , (c) t = 0.85 - 1 GeV 2 and (d) t = 1 - 2 GeV 2 . The parameters given in Table I are used. We 
compare in Fig. 5 (a) - (d) the calculated results for |-F(i)| 2 for the time-like momentum with the 
experimental data. It is necessary to take account of the to boson and p-Lu mixing to explain the dip 
in the form factor near the p meson mass. In Fig. 4 and Fig. 5 (a) - (d) the dashed curves denote 
the result obtained by taking account of p'" . We find that the assumption of p'" improves the result 
especially for the time-like momentum above 1.5 GeV 2 . The chi square is obtained as x 2 — 379.1 
for the case without p'" and x 2 = 355.1 for the case with p'". 



10 



4 Concluding remarks 



We analyzed the experimental data for the pion electromagnetic form factor by using the dispersion 
relation in which the QCD constraint is satisfied. The imaginary part of the form factor is given 
as the summation of the Breit-Wigncr formula with mixing among resonances and the QCD term. 
The latter is expressed as a power series in the running coupling constant multiplied by a function 
that assures the convergence of the superconvergence constraint. The function may be interpreted 
as representing the quark-photon and quark-hadron vertices. In order that the form factor satisfies 
the QCD prediction we imposed the asymptotic condition on the absorptive part, by which the 
coefficient of or in ImF® CD , that is cf D , is given in terms of the pion decay constant. The QCD 
constraint is very strong so that at least 0{a 2 R ) approximation is required to realize the experimental 
data. By taking the second order approximation for the QCD running coupling constant we are able 
to reproduce the experimental data of the pion electromagnetic form factor very well. The result is 
improved by taking account of the third order term, namely, c^ CD and c^ CD as free parameters. 
However, the experimental data are not sufficient to determine them accurately. 

We assumed the mixing among vector bosons. Especially mixing between uj and p(760) is nec- 
essary to realize the structure of form factor near the p meson mass as in Ref.fuJ. The structure 
is explained by the mixing very well as is illustrated in Fig. 5 (b). An enhancement is observed in 
the pion electromagnetic form factor for t = 1.1 ~ 1.3 GeV |l3| (see Fig. 5(d)), although there is no 
isovector vector boson with the mass between 1.2 and 1.7 GeV. To explain the enhancement, without 
assuming resonance with the mass about 1.2 GeV, large mixing between p' and p" is necessary. 

Analysis is performed by using the data of Amendolia et al.[[l4| and Bebek et al.[jl5| for the 
space-like momentum. As we have shown in Fig. 3 our result agrees with the recent experiments for 
the space- like momentum 

We examined the effect of unconfirmed vector boson p'" with the mass above 2 GeV. It is shown 
that p'" improves the result, especially for the time-like momentum with t above 1.5 GeV 2 (see Fig. 5 
(d)). 

The QCD scale parameter A is determined to be in the range 0.3 ~ 2 GeV from the experiments 
of the pion electromagnetic form factor; the best fit is attained for A ~ 0.8 GeV. A, obtained from 
the analysis of pion form factor, seems to be larger than that reported in the deep inelastic scattering. 
This implies that A becomes larger when the number of effective flavor is decreased [ Gl| . 
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Appendix A functional realization of theoretical results for 
the space-like momentum 



F f/1 = M " l + N lX + N 2 x 2 

U \t\ + M 2 1 + D x x + D 2 x 2 + D 3 x 3 ' 1 ' 



For the convenience of reproducing the theoretical form factor for the space-like momentum we 
express our calculated results by a functional approximation similar to the formula proposed by Gari 
and Kriimperman |TT 

' ,J 1 \ia; -t- iv 2 s 

l\X + D2X 2 + U 3 J 

where x = ln(l + \t\), with t being expressed in the unit of GeV 2 . M , a constant with the dimension 
of mass, is fixed at the p meson mass; M — 0.76 GeV. The parameters Ni and Di are determined 
so as to realize the calculated results given in Fig. 4. 

For the case without p" 1 , the parameters are determined as follows: N-y = 4.5919, N 2 = —0.62524, 
L>i = 4.9604, D 2 = -0.7049, D 3 = 0.0540. For the case with p" 1 the parameters are N x = 4.7597, 
N 2 = 0.30520, Di = 5.1071, D 2 = 0.3598, D 3 = 0.1904. The formula is valid up to \t\ < 10 GeV 2 ; 
the theoretical result is reproduced within the accuracy of 0.1% for the squared momentum \t\ < 10 
GeV 2 and for \t\ < 6 GeV 2 the error is less than 0.02%. 
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